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ABSTRACT 

The spectra of type Iln supernovae indicate the presence of apre-existing 
slow, dense circumstellar wind (CSW). If the CSW extends sufficiently far from 
the progenitor star, then dust formation shouldoccur in the wind. The light 
from the supernova explosion will scatter off this dust and produce a light 
echo. Continuum emission seen after the peak will have contributions from 
both this echo as well as from the shock of the ejecta colliding with the CSW, 
with a fundamental question of which source dominates the continuum. We 
calculate the brightness of the light echo as a function of time for a range of 
dust shell geometries, and use our calculations to fit to the light curves of SN 
1988Z and SN 1997ab, the two slowest declining Iln supernovae on record. We 
find that the light curves of both objects can be reproduced by the echo model. 
However, their rate of decay from peak, color at peak and their observed peak 
absolute magnitudes when considered together are inconsistent with the echo 
model. Furthermore, when the observed values of M B are corrected for the 
effects of dust scattering, the values obtained imply that these supernovae have 
unrealistically high luminosities. We conclude that light echoes cannot properly 
account for the slow decline seen in some Iln's, and that the shock interaction is 
likely to dominate the continuum emission. 
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1. Introduction 

Supernovae are designated to be of type Iln based on their spectra (Schelegel 1990). 
The group is quite heterogeneous, but a typical Iln would have the following properties: an 
Ha profile composed of a narrow peak (~ 100 km/s) sitting on a broad base (~ 10 4 km/s), 
no P Cygni absorption feature, a strong blue continuum, and slow spectral evolution. There 
is occasionally an intermediate width component (~ 2000 km/s) to the Ha line. Absorption 
lines tend to be weak or absent. The narrow Ha component is seen to vary in intensity 
and has thus been interpreted to indicate the presence of a pre-existing slow-moving stellar 
wind. These winds can be very dense, with M as high as 10~ 2 M /yr as is the case with 
SN1997ab (Salamanca et.al. 1998). The overall spectral properties of these supernovae 
have successfully been interpreted as being due to a shock interaction between the rapidly 
expanding supernova ejecta and this dense (10 6 — 10 8 cm~ 3 ) stellar wind (Chevalier & 
Fransson 1985, Chugai 1991, Terlevich 1994). Iln's comprise about 5% of observed type II 
supernovae. 

A number of Iln's show an abnormally slow decline in their light curves. SN1997ab 
and SN 1988Z are the extreme examples of this, both declining by only about 5 magnitudes 
in 1000 days. Simple models of the shock interaction suggest that the shock can produce 
sufficient optical continuum emission to account for this slow decline (Terlevich 1994, Plewa 
1995). It was suggested by Chugai(1992) that light echoes might also provide a significant 
contribution to the late time light curves of Iln's, i.e. the observed continuum should be a 
combination of shock and echo contributions. When the supernova explodes, it vaporizes 
any dust in the CSW out to a radius of ~ 10 16 cm (Pearce & Mayes 1986). Any dust that 
lies outside this evaporation radius will scatter the supernova light. Scattered light travels 
a longer path before reaching us and so is delayed relative to any unscattered light, i.e. a 
light echo is produced. As little as 1O~ 6 M of dust is required to make the circumstellar 
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shell of matter optically thick to dust scattering. 

We here investigate whether light echoes could help account for the slow decline in the 
light curve seen in some Iln's. We have performed a detailed calculation of the light echo 
produced by a range of dust shell geometries. Our assumptions and method of calculation 
are described in the next section. Our results are presented in Section 3, and in Section 4 
we fit to the observed light curves of SN 1988Z and SN 1997ab. Section 5 is devoted to the 
discussion of our results. 

2. Calculating the Light Echo 

The light echo in the U, B, V, and R bands was computed by a Monte Carlo 
simulation. The calculation requires knowledge of the dust size distribution, the dust spatial 
distribution, the optical properties of the dust, and the underlying Type II supernova light 
curve. 

The dust was assumed to be composed of a mixture of silicon and graphite grains, in 
the ratio 3:1 by number. Both grain types were assumed to have a grain size distribution 
described by the MRN grain model (e.g. Evans 1994), i.e. the number density of 
grains varies as n(a) oc a -3 ' 5 where a is the grain diameter. Our grain size ranged from 
a-min = 0.005/im up to a max = 1.0/im. 

We assumed that the dust lies in a spherical shell of inner radius R m in and outer radius 
Rmax about the star. R m i n is the radius out to which the supernova evapourates all the 
dust in the CSW. R m in can be estimated from first principles and for a type II supernova 
has a value of a few times 10 16 cm (Pearce & Mayes 1986). We used this value as a guide, 
and took both R m in and R max to be free parameters. We assumed that the dust density in 
the shell varies as p(r) oc r~ 2 . The amount of dust in the shell is characterized by r, the R 
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band optical depth for a purely radial photon trajectory. The total dust mass is related to 
r by the following expression: 

M dust = — (af^ ax - 0^f n ) RmaxRmin {f Pgr + (l ~ f)PSi) X 

T 



faZin ^ aa 15 Si=afes,scat {fQi,gr + (1 



Si J 



7.0 x HT 7 r | 1 [ ^ 1 M & 



^o.oipcy vi-0pc / 

where / = 0.25, p gr = 2.5gcm -3 and psi = 2.35gcm -3 are the densities of graphite and 
silicon respectively, and Q a b s and Q sca t are the absorption and scattering efficiencies. 

The photons are either absorbed by the dust, or they scatter elastically. The cross 
sections for the dust-photon interactions are functions of both the photon wavelength and 
the dust grain size. We used the scattering and absorption cross sections, cr a b s (= 47ra 2 Q a 6s) 
and a SC at{= 47ra 2 Q scai ), calculated by Draine(1987). The angular distribution for the 
scattering is given by the Henyey-Greenstein function (Henyey & Greenstein 1941), 

1 - q 2 

F(9) = y - m (2) 

(l + # 2 -2£cos6>) 3/2 

where g is the degree of forward scattering and is a function of wavelength and grain size. 
Values for g were also taken from Draine(1987). 

We used the composite II-L and II-P light curves calculated by Doggett and 
Branch(1985) as input. They provide averaged light curves for the B and V bands. We 
took the U and B bands to be initially identical, and likewise the V and R. The input 
supernova was taken to have B — V = at peak. 

EDITOR: PLACE FIGURE |T| HERE. 



Our Monte Carlo code follows the photon trajectories in three dimensions. It keeps 
track of the time delay for each photon, as well as the photon weight. The weight was 
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decreased at each interaction by the absorption probability. The photon was also split 
at each interaction into a part which has no further interactions and thus escapes, and a 
second part which is forced to interact at least once more. This substantially improves the 
output statistics. Our code was tested against a number of analytical test cases, and also 
against the calculations of Chevalier(1986), who calculated light echoes in the limit of low 
optical depth. 

Our code produces the R band light curve as well as the V — R, B — V, and U — B 
colors as output. The output is dependent on the choice of R m in, Rmax, and r. 10 6 input 
photons were used to calculate each band for each choice of the model parameters. Our 
output does not include any contribution from line emission. 

3. Results 

Our echo model has three parameters: R m i n , R m ax and r. Figures 2,3, and 4 show the 
effect of varying each of these in turn. 

EDITOR: PLACE FIGURE @ HERE. 
EDITOR: PLACE FIGURE g HERE. 

EDITOR: PLACE FIGURE | HERE. 

We see from the figures that the echo is fairly insensitive to both R m i n and R ma x as long 
as the shell is thicker than a few hundredths of a parsec. As the shell becomes thinner, 
the light curves fall off faster, and more closely trace the intrinsic light curve. We can also 
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see that the presence of a small amount of dust has a pronounced effect on the light curve 
at late times. The optical depths in the U, B, V and R bands are always in the ratio 
1.56:1.35:1.20:1.00. The V — R, B — V and U — B colors remain nearly constant after about 
200 days for all the non-thin shell cases. For the the thin shells the colors tend to redden at 
late times. 

While there are variations amongst models, the following relations serve as useful 
approximations: 

(U - B) peak ~ 0.17r , (B - V) peak ~ (V - R) peak ~ 0.15r 
Amu ~ l-4r , Arris ~ l-3r , Amy ~ l.lr, Amjj ~ 0.95r 
/3ioo,f/ ~ 4.3 - 1.5r, (3 W0 ,b ~ 4.3 - 1.4r , /9i 00 ,y ~ 4.0 - 1.2r, /3i o,r ~ 4.0 - l.lr 

Am is the number of magnitudes by which the peak of the light curve has faded due to 
dust scattering compared to the peak in the composite type II light curve, and /3 W0 is the 
number of magnitudes by which the light curve drops from peak over the first 100 days 
after peak. 



4. Fits to SN 1988Z and SN 1997ab 

SN 1988Z and SN 1997ab are the two slowest declining Iln's on record. (SN 1995N has 
a similar decay rate, but data on this object is not yet available). The slow decline suggests 
that these objects are embedded in a dense CSW. This is bourne out by analysis of the 
spectra of these objects. Salamanca et.al.(1998) calculate that the electron density in the 
CSW of SN 1997ab has a value close to 10 8 cm" 3 . Chugai(1992) found a similar result for 
SN 1988Z. We have used our echo model to fit to the light curves of these two supernovae. 
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4.1. SN 1988Z 

SN 1988Z was discovered on 1988 December 12 in the galaxy MGC+03-28-022 
(Cappelloro & Turatto, 1988). The supernova was past peak at this time, but its high 
luminosity (mg = 16.75, Mg = -18.3 for H = 65 kms _1 Mpc _1 ) argues that the maximum 
occured not long before discovery. The supernova faded monotonically at an unusually slow 
rate. Over the first 137 days it faded by only 1.5 magnitudes in both the B and V bands, 
whereas the composite II-L light curve drops by some 5 magnitudes over the same period. 
The color evolution of SN 1988Z is very unusual (Turatto et.al. 1993). Over the first 100 
days B — V decreases from 0.4 to 0.1 mag. The supernova then becomes redder, until the 
color reaches a constant value of about 0.65 mag after 700 days. 

Figure 5 shows our fit to the B and V band light curves of SN 1988Z. The B band fit 
corresponds to r = 2.0 ±0.2, R m in = 0.10 ±0.02 pc, R ma x = 0.35 ±0.05 pc. This corresponds 
to a dust mass in the shell of (eq. (1)) (4.9 ± 1.1) x 1O _6 M . The narrow component of 
the Ho; line is unresolved, so if we adopt a value of 100 kms -1 for the wind velocity, and 
assume a constant mass loss rate, we infer that the progenitor star emitted material for 
about 3400 years before going supernova at a rate of (1.3 ± 0.3) x 10~ 6 flj^f) M yr _1 , 
where fdust is the fraction of the wind mass which has formed into dust. The fit to the B 
band determines the fit to the V band, i.e. the V band is not fit separately. In both cases, 
the fit is reasonably good up to about 800 days, after which the observed light curves are 
seen to flatten significantly with respect to the model predictions. 



EDITOR: PLACE FIGURE | HERE. 
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4.2. 



SN 1997ab 



SN 1997ab was discovered on 1996 April 11 in the galaxy HS 0948+2018 (Hagen & 
Reimers 1997). It had a B magnitude of 14.7 which corresponds to an absolute magnitude, 
Mb = —19.1 (Hq = 65kms _1 Mpc~ 1 ). The date of peak is uncertain, but the high luminosity 
at discovery argues that the peak had occured recently. The light curve of this object is 
unfortunately not well sampled, there being only three other B band observations and two 
V band observation since discovery (Hagen, Engels & Reimers 1997, Schaefer & Roscherr 
1998, 1999). 

Our model fits the B band light curve (figure 6) for r = 3.0 ± 0.3, R m i n = 0.10 ± 0.02 
pc, Rmax — 0.30 ± 0.05 pc. This implies a dust mass in the shell of (6.3 ± 1.4) x 1O _6 M . 
The wind velocity is measured to be 90 kms -1 , and so we infer that the progenitor 
emitted material for some 3000 years before going supernova at an average rate of 



approximately 10 _2 MQyr _1 from the analysis of the spectrum of this object. They find that 
the emitted material lies in a shell that extends no more than 0.05pc from the progenitor. 
This implies that the progenitor had a brief episode of strong mass loss shortly before going 
supernova. The material emited during this episode lies too close to the star to allow any 
dust formed to survive the supernova explosion. The dust which is responsible for the echo 
must have been emitted during an earlier, more sustained period of lower mass loss. 




Salamanca et.al. (1998) find a mass loss rate of 



EDITOR: PLACE FIGURE | HERE. 
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5. Echo Model versus Shock Model 

If the progenitors of Type II supernovae emit material for more than a few hundered 
years prior to going supernova, then they should be surrounded by a dusty CSW. All the 
dust that lies beyond the evapouration radius of the supernova will scatter the supernova 
light and produce a light echo. There is likely to be an echo component in the light curves 
of many Type II supernovae, the only quesiton is how significant a component it is. 

We have seen in the previous section that the echo model can provide reasonably good 
fits to the light curves of SN 1988Z and SN 1997ab. Chugai(1992), however, produced 
a similar quality fit to the light curve of SN 1988Z with a simple model of the shock 
continuum emission. Fits to the light curves alone are thus not sufficient to discriminate 
between the models. 

In figure 7 we have plotted(S — V) pea K against /9ioo,b- The region between the two 
solid lines is the region of parameter space allowed by the echo model. As the amount 
of circumstellar dust increases, we expect the supernova light to suffer more scatterings, 
producing an echo that becomes progressively redder and longer lived. In shock models the 
optical light is primarily produced by the reprocessing of X-rays. The color of the optical 
light produced in this way is not very sensitive to the amount of circumstellar material 
present, so on the plot of (B — V) pea k versus /3ioo,b ; we do not expect to see a correlation. 
SN 1988Z and SN 1997ab both lie within the region allowed by the echo model. The other 
Iln's plotted on the figure, however, do not. Their distribution shows no definite correlation, 
and is thus consistent with the shock model prediction. The Iln's plotted represent the 
complete sample for which (B — V) pea k, Pwo,b, and MB, pea k are available in the literature 
(Patat et.al. 1994). The data for SN 1998S was taken from the CfA Supernova Group 
website. 

In figure 8 we have plotted the observed Ms tP eak versus Pwo,b- In the echo model, for a 
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given underlying supernova luminosity, as the dust mass increases, the decay rate decreases 
and the observed M B ^ peak increases, i.e. the supernova appears fainter as more of the 
emission at peak is scattered by the dust. The opposite trend is expected for shock models, 
i.e. we expect the peak luminosity to increase as the amount of circumstellar material 
increases. In an M B ^ peak versus /?ioo,b plot, shock models thus predict a downward sloping 
line. The observed Iln points are consistent with the shock prediction, and inconsistent 
with the prediction of the echo model. In particular, the points for both SN 1997ab and SN 
1988Z lie outside the region allowed by the echo model. 

In figure 9 we have plotted the observed MB, pea k versus (B — V) pea k- The echo model 
predicts that the emission should redden as the amount of dust, and thus Mb, peak-, increases. 
The observed data points are inconsistent with this prediction. 

Another argument against the echo model is the implied values for MB, pea k for the 
underlying supernovae. SN 1988Z had an observed M B - pe ak of -18.3 and SN 1997ab, -19.1. 
If we use the values of r, R m %n and R ma x obtained from the fits to the light curves to find 
Arris, we find that the underlying supernovae must have had MB tPea k of -20.9 in the case of 
SN 1988Z, and -22.8 for SN 1997ab. This would make these supernovae by far the brightest 
Type II supernovae yet observed. 

We conclude that light echoes cannot properly account for the slow decline seen in some 
Iln's. The predictions of shock models are consistent with the data, and shocks are likely 
to dominate the continuum emission. The absence of a strong echo component in the light 
curves of Iln's argues that the progenitors of these supernovae undergo the bulk of their 
mass loss just prior to going supernova, and thus very little dust survives the explosion. 
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Fig. I. — (a) and (b) respectively show the composite II-L and the composite II-P light 
curves used as input for the echo calculations, (c) and (d) show the values of B — V for each 
set of composite curves Doggert & Branch(1985). 

Fig. 2. — The effect on the light curves due to changes in r. R m i n and R max are held fixed 
at 0.1 and 1.0 pc respectively, (a) and (b) show the R band light curve for input II-L and 
II-P respectively. The remaining plots are for V — R, B — V and U — B for each of the two 
sets of input. The curves have been shifted downward by A mags for clarity. 

Fig. 3. — The effect of changes in R min . The curves are for R min = 0.01, 0.5, 0.95 and 0.99 
pc. r is held fixed at 1.0 and R max a t 1-0 pc. (a) and (b) show the R band light curve for 
input II-L and II-P respectively. The remaining plots are for V — R, B — V and U — B for 
each of the two sets of input. The curves have been shifted downward by A mags for clarity. 

Fig. 4. — The effect of changes in R ma x- The curves are for R ma x = 0.11, 0.15, 0.5, 1.0 and 
2.0 pc. r is held fixed at 1.0 and R m in at 0.1 pc. (a) and (b) show the R band light curve 
for input II-L and II-P respectively. The remaining plots are for V — R, B — V and U — B 
for each of the two sets of input. The curves have been shifted downward by A mags for 
clarity. 

Fig. 5. — Fit to the B (crosses) and V (circles) band light curves of SN 1988Z (Turatto et.al. 
1993). 

Fig. 6. — Fit to the B (crosses) and V (circles) band light curves of SN 1997ab. 

Fig. 7. — Plot of (B — V) pea k versus /?ioo,b- The area between the solid lines are the values 
allowed by the echo model. The underlying supernova was taken to have B — V — at peak. 

Fig. 8. — Plot of M B>peak versus /?ioo,b- M BtPeak for the underlying supernova was allowed to 
range between -16 and -19. 
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Fig. 9.— Plot of M B:Peak versus (B - V) peak . For the underlying supernova, M B:Peak was 
allowed to range between -16 and -19, and we assumed (B — V) peak = 0. 
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